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Abstract. In this note, we study energy estimates for Einstein vacuum equations in 
order to prove the formation of black holes along evolutions. The novelty of the paper 
is that, we completely avoid using rotation vector fields to establish the global exis- 
tence theorem of the solution. More precisely, we use only canonical null directions as 
commutators to derive energy estimates at the level of one derivatives of null curvature 
components. We show that, thanks to the special cancelations coming from the null 
structure of non-linear interactions, desirable estimates on curvatures can be derived 
under the short pulse ansatz due to Klainerman and Rodnianski [3] (which is originally 
discovered by Christodoulou [2 ). 



1.1. A Brief History. Penrose singularity tlieorem states tliat if in addition to the domi- 
nant energy condition, the space-time has a trapped surface, then the space-time contains 
singularities. The weak cosmic censorship conjecture asserts that under reasonable phys- 
ical assumptions, singularities should be hidden from an observer at infinity by the event 
horizon of a black hole. Thus, by combining these two claims, to predict the existence of 
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black holes, it suffices to exhibit one trapped surface in a space-time,. In other words, al- 
though many supplementary conditions are required, we regard the existence of a trapped 
surface as the presence of a black hole. 

A major challenge in general relativity is to understand how trapped surfaces actually 
form due to the focusing of gravitational waves. In a recent breakthrough [2], Demetrios 
Christodoulou gave an answer to this long standing problem. He discovered a mechanism 
which is responsible for the dynamical formation of trapped surfaces in vacuum space- 
times. In the monograph [2|, in addition to the Minkowskian fiat data on a incoming 
null hypersurface, Christodoulou identified an open set of initial data (this is the short 
pulse ansatz) on a outgoing null hypersurfaces. Based on the techniques developed by 
himself and Klainerman in the proof of the global stability of the Minkowski space-times 
[3], he managed to understand the whole picture of how various geometric quantities 
interact along the evolution. Once the estimates on curvatures are established in a large 
region of the space-time, the actual formation of trapped surfaces is easy to demonstrate. 
Christodoulou also proved a version of the same result for the short pulse data prescribed 
on past null infinity. This miraculous work provides the ffist global large data result 
in general relativity (without symmetry assumptions) and opens the gate for many new 
developments on dynamical problems related to black holes. 

In Klainerman and Rodnianski extended aforementioned result of Christodoulou. 
They significantly simplified the proof of Christodoulou (from about six hundred pages 
to one hundred and twenty). They also enlarged the admissible set of initial conditions 
and show that the corresponding propagation estimates of connection coefficients and 
curvatures are much easier to derive. The relaxation of the propagation estimates are just 
enough to guarantee that a trapped surface still forms. Based on the trace estimates devel- 
oped in a sequence of work of the authors towards the critical local well-posedness for Ein- 
stein vacuum equations, they reduced the number of derivatives needed of Christodoulou 
in the argument from two derivatives of the curvature (in Christodoulou's proof) to just 
one. More importantly, Klainerman and Rodnianski introduced a parabolic scaling in |1] 
which is incorporated into Lebesgue norms and Sobolev norms. These new techniques 
allow them to capture the hidden smallness of the nonlinear interactions among different 
small or large components of various geometric objects. The result of Klainerman and 
Rodnianski can be easily localized with respect to angular sectors, has the potential for 
further developments, see [3]. We remark that Klainerman and Rodnianski only concen- 
trated on the problem on a finite region. The question from past null infinity can be 
solved in a similar manner as in [2] once one understand the picture on a finite region. 
The problem from past null infinity has been studied in a recent work by Reiterer and 
Trubowitz, [6]. 

1.2. Novelty of the Paper. One common feature of the proofs in [2] and [1] is that, in 
order to derive energy estimates on one or higher derivatives of curvature components, 
they all constructed three angular momentum vector fields 0^^\ O^"^^ and O*-^^ which 
essentially captured almost rotational symmetry of the space-time. As far as the author 
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aware, the main reason of using 0^*^'s is that the energy estimates on {Co(i)R)ap-iS behave 
well because one can take advantage of cancelations from the pseudo-symmetry of 0^*^'s. 
Here, the £o(0 is the modified Lie derivative defined in [3] and Rap-yS is the curvature of 
the space-time. We can observe this advantage in the proof of stability of Minknowski 
space-time |3] where (Co(i)R)ai3'y& yields better decay estimates. 

There is one obvious defect of the modified Lie derivative. As usual Lie derivatives, the 
Co(i) is not tensorial in O*-*-*. In fact, it evolves one derivative of 0^*\ In other words, if 
we use modified Lie derivatives, we may lose immediately one derivative. We have two 
remedies to this loss of derivatives: in [2], one relies on higher order derivative estimates; 
in [1], one makes use of more subtle trace estimates. 

The above discussion can be summarized as follows: roughly speaking, a good estimate 
on Lie derivatives relies on the almost symmetries 0*^*^'s, but Lie derivatives causes a loss 
of derivative as we just explained. Hence, to avoid this loss, we shall give up the use of 
the pseudo-symmetries. 

In this paper, we propose an approach to derive energy estimates on curvatures without 
constructing rotational vector fields at all. In particular, in stead of using modified Lie 
derivatives, we work with covariant derivatives to save one derivative. Of course, we have 
to pay a a price of controlling much more error terms. Nevertheless, this significantly 
simplifies the proofs compared to either [2] or [1]. In particular, compared to [2], we use 
only one derivative in curvature; compared to |1] , instead of using trace inequality, all the 
estimates are derived from the classical Sobolev inequalities. 

We also want to mention that in the thesis of L. Bieri, see [1], based on a more general 
asymptotic assumptions, she gave a simplified proof of the stability of Minkowski space- 
time. She managed to derive decay from the time vector field and the conformal scaling of 
the space-time which allowed her to circumvent rotational vector fields in that situation. 
But the current situation is different from [1| : we do not use Lie derivatives at all and the 
lower regularity forces us to explore more structures from the Einstein equations. 

1.3. Structure of the Proof. The main observation arises from the second Bianchi 
identities. Roughly speaking, they explicitly show how one expresses angular or rotational 
derivatives of curvature components V^E' in terms of some null derivatives of curvature 
components plus lower order non-linear terms. Namely, they can be written schematically 
as 

Vl^ = V^ + /.o.t., 

or 

Vl^ = V^ + /.o.t., 

where we use lower order terms l.o.t. to collect all nonlinear interactions and L, L are two 
standard null directions under the framework of double null foliations. Thus, up to lower 
order corrections, to obtain the estimates on V\E' rotational derivatives on curvatures, it 
suffices to control Vl\E' or V^^^ null derivatives of curvatures. Thus, we identify our main 
targets to be D^Ra^-ys and Dj^Rap-ys and we shall derive energy estimates for them. 
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According the above idea, we use the modified short pulse ansatz proposed Klainerman 
and Rodnianski in We remark that this ansatz allows more large components than the 
original short pulse data discovered by Christodoulou in |2]. We shall derive the energy 
estimates based on Bel- Robinson tensors associated to DiRajB-ys and DiRai^^s- To deal 
with error terms, namely terms J, J and K in the proof in following sections, we have 
to take account of the special structure of those terms. In reality, some generic terms in 
error term may cause a loss of 6^2 which prevent us from closing the bootstrap argument. 
To avoid this loss, there are typically three techniques to use: 

1) We bound a product of two term in L^^^^ by two L^sc) estimates on each term 
instead of one L^^-j estimate on one of them and one L^^^^ estimate on the other. 

In most of the situation, this trick saves a 6^. 

2) When we integrate a product of terms on some null hypersurface or on a domain 
in the space-time, we use integration by parts to move a bad derivative, typically 
V3 or V4, from one term (for whom this bad derivative may cause a loss of 5~2) 
to another (for whom this bad derivative is not bad at all, namely, there is no loss 
in S). Combined with Bianchi identities, this procedure may save a 5 2. 

3) For some generic term in error estimates mentioned above, though it appears that 
it causes a loss of S (which can not be retrieved by using the trick 1) and 2)), we 
can use in fact either signature considerations or a precise computations to show 
that this term does not show up at all in the error estimates. This manifests the 
special cancelations in the error termsj^ 

The whole proof is to combine these three tricks. The paper is organized as follows: 
in next section, we recall basic definitions and estimates from ^ and we state the main 
theorem; in following sections, we derive energy estimates for derivatives of curvature 
components in the following order: V4a, Vsa and then V4\E'i for ^ a and Va^E' for 
^ ^ a. 



2. Main Result 

2.1. The Double Null Foliation Framework. We briefiy recall the double null foli- 
ation formalism. See [2] for more precise definitions. We use V = S) to denote the 
underlying space-time and use g to denote the background metric. We assume that V is 
spanned by a double null foliation generated by two optical functions u and u and they 
increase towards the future, < u < 1 and < u < 6. We use Hu / to denote the 
outgoing / incoming null hypersurfaces generated by the level surfaces of u / u. We use 
Su,u to denote the space-like two surface Hu H We denote by H^r^'-^^ the region of 
Hu defined hj Ui < u< U2'i similarly, we can define Hj;^^''^^\ 



This cancelation can also be observed much more directly by another way of deriving energy 
estimates, namely, multiplying Bianchi identities and integrating directly on a give domain. The author 
would like to thank Igor Rodnianski for communicating this idea. 
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The shaded region on the right represents the 
domain V{u,6) with < m < 1. The func- 
tion u is in fact defined from —1 to 5. When 
u < 0, this part of Hq is assumed to be a flat 
hght cone in Minkowski space with vertex lo- 
cated at u = —1. We use ro ~ 2 to measure 
the maximal radius of the flat part of Hq. In 
the trapped surface forms at m ~ 1 and 
u = S. 



Let (L, L) be the null geodesic generators of the double null foliation and we define 

2 

the lapse function Q by g{L,L) = —7^77. The normalized null pair (63,64) is defined by 

63 = fiL, 64 = fiL, ^'(63,64) = —2. On sphere Suu we choose an arbitrary orthonormal 
frame (61,62). We call 

We use D to denote Levi-Civita connection defined by the metric g and we define the 
connection coefficients as follows. 



X.ab 



g{Dbe4,eb), Va 



^g{D3ea,ei), w = -^^(D463, 64), 



g{Dbe3,eb), = --giD^Ca^Cs) 



1 1 

^= -^fi'(^3e4,e3), (a = ^g{Daei,e3), 



where Da = De^. On Su.u, V is the induced connection; V3 and V4 are the projections 
to Su.u of the covariant derivatives D3 and D4. 

Given a Weyl field W, we introduce its null decomposition with respect to the given 
null frame, 



a(W^)ab = W{ea, 64, 6fc, 64), f3{W)a 
a{W)ab = W{ea, 63, 6fe, 63) (3iW)a = ^W{ea, 63, 63, 64 



^W{ea, 64, 63, 64), p{W) = |IV(64, 63, 64, 63), 



a{W) = -*W{e,,e^,e,,es) 



where *Vr is the space-time Hodge dual of W. When W is the Weyl curvature tensor, we 
use a, a, (3, (3, p, a to denote its null components. 

We recall the null structure equations for the Einstein vacuum space-times (see [4]). 
The originally Einstein field equations are 



0, 



We use Greek letters a, /3, • • • to denote an index from 1 to 4 and Latin letters a, 6, • • • to denote 
an index from 1 to 2. We also N to denote the null direction either L or L. Repeated indices are always 
understood as taking sums 
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where Rap is the Ricci curvature of the underlying space-times. We express this tensorial 
equations by using the null frame, by definition, this yields the null structure equations. 
We only list the transport type null structure equations which are relevant to the current 
work. 



V4trx + ^(trx)^ = - \x\^ - 2utTX, V4X + trxx = -2wX - «, 

Vstrx + ^(trx)^ = - \x\^ - 2urtYx, V3X + tr^X = -2wx - a, 

WaT] = -X- {t] - Ti) - (3, V3?7 = -X- iv- V) + 

3 1 11 

V4W = 2ww + -|r7 -r/p - -iv-v) ' iv + v) - gl^ + !Zl^ + 3^' 

3 1 11 

V3W = 2uuj + -\v-v\^ + ^iv-v)-iv + v)- ^\v + !l\'^+ 

V4trx+ ^trxtrx = 2wtrx + 2divr7 + 2|r7p + 2p - x ■ X, 
Vatrx + -trxtrx = 2a;trx + 2divr7 + 2\ri\'^ + 2p — x ■ X 
V4X + ^trxx = V®?7 + 2ux - ^tixx + V^V, 
V3X + ^trXX = V§?7 + 2u;x - ^trxx + V^V- 



2.1 

2.2 
2.3 

2.4 
2.5 
2.6 
2.7 
2.8 
2.9 



We also express the Bianchi equations relative to the null frame to derive null Bianchi 
equations]^ 

V3a + ^trxa = V§/3 + 4ua - 3{xp +*xa) + (C + 4r7)§/3, (2.10 

2.11 
2.12 



V4/3 + 2trx/3 = div a - 2w/3 + r] ■ a, 
V3/3 + trx/3 = Vp +*Va + 2w/3 + 2x-/3 + 3(r/p +*r]a), 

+ hixa = -div*/3 + 1^ .*« 277 ■*/?, 

V3a + ^trxa = -div *P + ^x *(l -C-*^- 2r] ■* (3_, 
V4P + ^trxp = div/3 - ■ « + C ■ /3 + 2r7 ■ /3, 

V3P + ^trxp = -div/3 - ^x ■ a + C ■ /3 - 2r/ ■ /3, 

V4^ + trx^ = -Vp +* Va + 2u;/3 + 2x ■ /3 - 3(r7p -*r/cT), 
V3/3 + 2trx/3 = -diva - 2ujj3 + r] ■ a, 



2.13 
2.14 
2.15 

2.16 

2.17 
2.18 



V^a + ^trxa = -Vg/3 + Aua - 3(xp -*X(^) + (C - 4r7)®/3, (2.19) 



div is the divergence on Su,u and Div is the space-time divergence. 
We can eUminate C by C = I (*? — ^) ■ 
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2.2. Energy Estimates Scheme. We review our scheme for energy estimates on Weyl 
fields, see [3] for the original resource. Assume a Weyl field Wa/S'yS and its Hodge dual 
*Wa/3'y5 solve following divergence equations with source terms 

D^Wa/B-fS = Jl3^S, D°'*Wal3-fS = *Jp'r5, (2.20) 

Jai3'y and *Ji3-ys are called Weyl currents. 

Remark 2.1. For vacuum, the curvature tensor Ra/s-yS is a Weyl field with zero currents 

D'^Rap-yS = 0, D°'*Rai3-yS = 0. (2-21) 

The Bel- Robinson tensor associated to Wap-yS is defined as follows]^ 

It is fully symmetric and traceless in all pair of indices. Moreover, it satisfies the dominant 
energy condition which allows one to recover estimates for Weyl field W. In pratical terms, 
this condition can be expressed by formulas, 

Q4444 = 2|«p, Q3333 = 2|a|^ Q4443 = 4|/3|2, Q3334 = 4|/3|2, Q4433 = 4(p2 + fx^). 

We also list other null components of Q for future use, 

Qa444 = ^OlabPbi Qa333 = ~"4a^^^^, 

Qa344 = 4p/3, - 4a Qa433 = '^P^^ ' 4cr (2.22) 
Qafe44 = 2|/3|2 + 2pa - 2a *a, Qab^z = 2|/3p + 2p« + 2a*«, 
Qabsi = + 2(p2 + a^)5ab- 



In view of (2.20), Q enjoys the following divergence equations 

= ^/3"/^^7- + Wp^.'^J^s. + W/^'^^^, + W//*J^5.. (2.23) 

Given vector fields X, Y and Z, we define the current P\W]{X.,Y., Z)^ associated to X, 
y, Z and W by 

= P[W]{X, r, Z)^ = Q^p^sX^Y^'Z'. 

Thus, 

D^P^ = D^Q^xYZ + (vr ■ Q){X, Y, Z), (2.24) 
where ^^^tt is the deformation tensor of X defined by ^^''rcap = ^{DaXj^ + Di^Xa) and 

(vr ■ g)(X, Y, Z) = Q^p./''\^^Y^Z' + Q^^^^^^V^^Z^X^ + Q^p,s^^\^P X-^Y' . 

We integrate (2.24) on V{u,u) to derive the fundamental energy identity 



^ We shall use short hand notations Q for a for a{W), jS for /3(VF), p for p(Vt^), if there is 

no confusion in the context. 

^ L and L are corresponding normals of the null hypersurfaces and H_^. 
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I Q{X,Y,Z,L)+ I Q{X,Y,Z,L) 
^ \ =1 Q{X,Y,Z,L)+ f QiX,Y,Z,L) (2.25) 
4-// DiYQ{X,Y,Z)+ [[ {n-Q){X,Y,Z). 

JJv(u,u) JJv(u'a) 



'>{u,u) JJT>{u^u) 

We list non-zero components of deformation tensors of L and L as well as the non-zero 
component of D^I/ and D^I/ for future us^ 

(^Vgs = -8^]-V, ^"-W = 2^]-lr/^ (^V,, = 

= 2a;, D^L'^ = = -^-^7]", D^L^ = 

We consider one derivative of curvature Rap-yS in a null direction A^. One easy but 
important observations is that D^Ra^^s is still a Weyl field. We can commute with 
(|2^ to derivfl 

D'^DNRaji-yS = J p^Si D'^*DNRap-yS = *J li^Si (2.26) 

where 

^/375 ~ R^ N Rnu-fS + R^N^'^RfiPuS + R^ Ns'^ Rfi/B-yu + D^R^p^s-, (2.27) 

and 

= R^ N p'*Riiv'yS + R^ N-i^*RiiPv5 + R^ Ns'*Riip-iu + D^N'^D^*Rfj_j3^s. (2.28) 

2.3. Short Pulse Ansatz and Scale Invariant Formulation. We briefly recall the 
notions of signature and sca/e introduced by Klainerman and Rodnianski in [1]. Let 
be either a null component of curvature or a connection coefficient, we use Na{(l)), A'3(0) 
and N4^{4>) to denote the number of times (ea)j=i,2, respectively 63 and 64 appearing in 
the definition of 0. The signature of 0, sgn^cj)), and the sca/e of 0, sc(0), are defined as 

sgn{^) = 1 X Ar4(0) + ^ x Ar,(0) + x Ar3(0) - 1, sc(0) = -sgn{^) + ^. 

We list the signatures and scales for all connection coefficients and curvature components, 

''' We do not use 63 and 64 in order to avoid V4W and V3W which do not have certain estimates, 
see H] for details. 

® Recall that Dpf commutes with Hodge * operator. 
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3 
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2 
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2 
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3 
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-1 
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2 









We impose 



1 



sgn{Vi(t)) = sgn{(f)) + 1, sgn{V(j)) = sgn{(t)) + -, sgn{V3(f)) = sgn{(f)) + 0, 

sgn{(j)i ■ 02) = sgn{(t)i) + sgn{(t)2). 

We define the scale invariant norms for (p. Along null hyp ersurf aces Hu'^^ or H_l^'^\ 
On a two dimensional surface Su,u, 

^-sc{</.)-i 



Those norms are obviously related by formulas 



5 



-1 



(sc) ^ '— ' 



Those scale invariant norms come up naturally with a small parameter 6. Roughly speak- 
ing, it captures the smallness of the non-linear interaction. We have Holder's inequality 
in scale invariant form, 



1101 ■ 02||l^^^j(5„,^) < s'^llMill^is^JlMLll^is^,^) with - = - + 



1 _ 1 1 

P Pi P2 



(2.29) 



Similar estimates hold along null hypersurfaces. 

Remark 2.2. The rule of thumb for treating the nonlinear terms is, whenever one has a 



product of two terms, (2.29) gains a 55. We do have cases that (2.29) does not gain any 
power in 5. In fact, if f is a hounded (in usual sense) scalar function (say, hounded hy a 
universal constant), the hest we can hope is ||/-0||lp ^{Suu) ~ II0II-L|' ^{Suu)- -^'^ particular, 

^ - 1, we have to pay 



in this paper, for f can he trx = trx + trx^ where trx^ = 2ro+u-u 
special attentions to the appearance of trx, see [4j for more detailed descriptions. 
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We introduce a family of scale invariant norms for connection coefficients where p = 2, 4 
or oo, 1^ 

1^) — II (x? x) II Lf , + II (ti'X? '-^5 ^5 ^7 trx, cu) ll/^p 

{sc) (sc; 

Oi,p{u,u) = ||V(x,trx,w,r/,?7,x,trx,w)||iP^^(„_^),p ^ oo, 

= ||V^(x,trx,w,r7,r7,x,trx,w)|| 2 ,^{o.u)., 
as well as for curvature components, 

TZo{u,u) = 6^\\a\\ 2 .„(o,ii). + ||(/3,p,a,/3)|| 2 ,rr{o,u)., 

(sc) ^ ' (sc) V / 

T^oiu^U) = d^WPW ^ ,j:r(0,«)-, + ||(p,Cr,/3,a)|| 2 , (0,.) , 

7^i(u,m) = 55||V4a|| 2 .^(0,3i)s + ||V(a,/3,P, Cr,/3)|| 2 /rrCCM)-,, 

(sc) V ^ / (sc) V / 

= 5^||V3«|| 2 .r.(0,^)s + ||V(/3,P,Cr,/3,a)|| 2 /r^CO,-)-,- 

Finally, we introduce total norms. We define Co,4 = sup„ „ Oo,4(m, m) and T^o = 
sup„ „ 7^o(w, m); similarly, we can define supremum norms for other scale invariant norms. 
The total norms are defined as follows, 

7^ = 7^o + 7^l, 1i = 11q + 11i. 

We use and 7^(°)to denote total norms on the initial hypersurface Hq. 

In the above definitions, all norms are scale invariant except for ||x||/^2 (uu)^ IIxIIl^ (uu)^ 
(h )i llVsalUa (H ) and 11/3 La ,^ Those terms are understood to cause a loss 
of 5~2. Notice also /3 on incoming hypersurfaces is scale invariant. By abuse of lan- 
guage, we call those terms anomalies or anomalous if they cause a loss of 5~p in L^^^^ 
norm. Notice also all the connection coefficients are not anomalous in L^^^ norms. In- 
spired by this, we use ipg ('g' for good) to denote some (good) connection coefficient in 
{trx, cj, 7], ?7, trx, u} and 4) to denote an arbitrary connection coefficient. We use '^g to de- 
note some curvature component in {/?, p, cr, /3, a} and ^ to denote an arbitrary connection 
coefficient. We shall also put a 'g' as a lower index to other quantities in order to indicate 
that this quantity is not anomalous. For example, we can write '^{DaR)g = a{DaR) 
because \\a{DaR)h2 (H^o,n). < 1. 



9 We use shorthand notations IK?/-, V-', V'", ' ' ' )ll = M + IIV-'ll + WW + 
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2.4. Main Result. For Einstein equations witli cliaracteristic data prescribed on H^q 
(where tlie data is trivial) and Hq, we can freely specify the conformal geometry on Hq. 
In other words, we can specify x freely along Hq to fix an initial data for the evolutionary 
problem. We remark that, in contrast to the case where the initial data is give on a 
space-like hypersurface, there is no constraints and the data can be freely specified. 
We require the initial data x subject to the following norm is finite, 

2 

Io = S'^\\x\\l^{Ho)+ sup [X1'^^II(^^4)^X||l2(5o.s) 



0<u<5 



k=0 



A:=0 m=0 

i.e. our main assumption on the initial data is the following short pulse ansatz, 

Jo < oo. (2.30) 

Notice that the size of x can be as large as 5^^ so the problem is far away from small 
data regime. 



The ansatz (2.30) was introduced by Klainerman and Rodnianski in This initial 

data set is larger than those of Christodoulou's original short pulses. In fact, (|2.30) allows 

1 

more components than those of Christodoulou's to be as large as 6~2. In this ansatz, we 



allow V behaves as 6' 



; in Christodoulou's ansatz, V behaves as 1. 



Under the ansatz (2.30), with the help of null structure equations, we can easily derive 
the following estimates on initial outgoing surfaces Hq, 



Lemma 2.3. Under the ansatz (2.30), along the initial outgoing hypersurface Hq, if 5 is 



sufficiently small, there is a constant c{Iq) depending only on Iq, such that 

0(0) + 7^(0) <^(j^)_ 



(2.31) 



Thanks to this proposition, we shall replace (2.30) by (2.31). And we omit the proof 



and refer the reader to or Chapter 2 of [2J. The next proposition claims that we 
can control connection coefficients provided that we have bound on curvatures. This is 
Theorem A in ^ and we omit the proof. 

Proposition 2.4. Assume that TZ and TZ are finite in T>{1,S). Then there is a 

constant C depending only on 0^^\TZ and TZ such that^^ 

< C. (2.32) 
We now state our main theorem. This is the following propagation estimates, which 



asserts that if (2.31) holds on initial hypersurface Hq, thus on the whole T>{1,5) we can 



bound the curvature norms TZ and ^ by a function depending only on initial data. 



10 



Throughout the paper, we use C to denote a constant depending only on 0^^\TZ and TZ. 
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Main Estimates. Assume the short pulse ansatz (2.30) hence (2.31), thus if 5 is suffi- 
ciently small, we have 

n + ]Z<c{Io), (2.33) 
where c(Xo) is a constant depending only on the size of the initial data. 

The main consequence of our estimates is the following global existence theorem, 



Existence Theorem. Given initial data x satisfying (2.30), if S is sufficiently small, we 
can construct an unique solution of the Einstein vacuum equations 

Rap = 0, 

on the whole region V{1,6). 

Once we have the Main Estimates, the proof of the theorem is a typical continuity 
argument. We refer the reader to Chapter 16 of [2] for the detailed proof. We shall not 
repeat this argument. 

We want to remark that the formation of trapped surfaces can not be derived from the 
initial data ansatz (2.30). We have to add two more modifications to make sure that Hq 
is free of trapped surface and S{1, 6) is a trapped surface. More precisely, in addition to 
(2.31), we also assume that 



sup X^5^||(5^)'V'xI|l2(5o,j 



0<u<5 



k=2 



for sufficiently small e such that < 6 <^ e and we also assume x satisfies. 



2(ro 



u 



< 



\x{0,u)\'^du< 



2(ro - 5) 



(2.34) 



(2.35) 



where Cq is a universal constant and ro ~ 2 measure the maximal radius of the fiat part 
of Hq. These condition guarantees the dynamical formation of trapped surfaces, we refer 
the readers to [1] for details. 

2.5. The Bootstrap Argument. We use a bootstrap argument to prove the Main 
Estimates. To be more precisely, we assume that TZ and TZ are finite. This assumption 
holds initially near Hq. Our main task is to carry out the following estimates under this 
assumption, 

n + n<c{io){i + {n + nfs) + c6^. (2.36) 

Thus, if 6 is sufficiently small, we have 

n + n<c{io). 



This yields the Main Estimates ( |2.33[ ). 

The derivation for (2.36) consists of two steps. The first step is to derive energy es- 
timates for curvature components; the second step is to derive energy estimates for one 
derivatives of curvature components. The first step is done in [Ij: 
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Proposition 2.5 (Proposition 14.9 in [4j). If 6 is sufficiently small, we have 

||(/3,p,a,/3)|L2 (^^) + ||(p,a,/3,a)|L2 < Xq + c(Xo)7^^ + C^i. 
As we explained in the introduction, our main target is the second step, i.e. to derive 



energy estimates for one derivatives of curvature components in the form of (2.36) without 



constructing rotational vector fields 0''*'''s. This is done in the rest of the paper. 

3. Preliminary Estimates 

In this section, we collect some estimates either already derived in ^ or relatively easy 
to prove directly. 

3.1. Improved Estimates on Curvature. We have estimates for curvature on each 
leaf S = Su^u- 

Proposition 3.1 (Proposition 6.6, Proposition 6.9 in // 5 is sufficiently small, we 
have 

,(s) + ll(/3,P,c^,/3,a)ilL2 AS)^C. 

(sc) ^ — (sc) ^ ' 

We also have L^^^^ estimates on curvature components. We recall again that the con- 
stant C depends on TZ and TZ which has information on the one horizontal derivatives of 
the curvature. 

Proposition 3.2 (Lemma 10.1 in fl]). If 5 is sufficiently small, we have 

^H(^\\Lj AS) + ll(/3,P,fx,/3,a)||i4 IS) ^ C. 

ysc) l^cj 



This proposition is proved by virtue of null Bianchi equations (2. 10)- (2. 19) and following 



scale invariant Sobolev trace type inequalities stated as Proposition 4.15 in ^4]. 
Proposition 3.3. Given an arbitrary tensor field (p, we have 

llf^^^iS) < i^^'Uhl^^iH^) + \\V<P\\Ll^^^H^))ks'^U\\Ll^^iH^) + \\^4<P\\Ll^^iH^))K 
L^^^iS) < (S'^UWLl^.iHJ + I|V0||l2^^^(Hj)^(5^||0||l2^^,(HJ + W^S^Ll^.iHjy^- 

Integrating along null hypersurfaces, we derive 
Corollary 3.4. If 6 is sufficiently small, we have 

+ II(/3,P,c^,/3,«)IIl4^^)(h„) < C, 
^''h\\Ll^,{Hj + ||(/3,p,a,^,a)|U4^^^(^,j < C. 
We also collect some more precise estimates on x, x and u. They are essentially 



developed in Section 5 of |1] provided we bound a by Proposition 2.5 
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Proposition 3.5. If 6 is sufficiently small, we have 



3.2. Comparison Estimates. The aim is to compare null components of derivatives of 
curvature and derivatives of null components of curvature. The first lemma compares 
VAr\l' with "^{DnR) where N = L or L. 

Lemma 3.6. For N = L or L, we have 

Proof. By direct computations, in schematic form, we have the following commutator 
formula]^ 

^pTvi?) - ViV^ = ■ ^ 
To estimate V'g ■ ^ on some null hypersurface (say H without loss of generality), we use 
Ll estimates on both factors and notice that \E' may cause a loss of 5*, thus we have 



since llV^gllL* (m ~ C and ||\I'|U4 (jj) < C5^^. This completes the proof. □ 

The next lemma, by investigating null Bianchi equations, compares VAr\E' with V\i/'. 

Lemma 3.7. We have following estimates, 

(1) . ForV^Oi, we have 

(sc) ^ (sc) ^ — Vl' 

(2) . For V 4- {Vsa, V4a, Vsa}, except V4/3 on if^ anc? V3/3 on H^, we have 

(sc) ^ (sc) ^ — 1^' 

Remark 3.8. We can not estimate on ||V4/3||j;^2 ^ or ||V3/3||£2 (j:^^-, ^/ie moment. 
The reason is that we do not have any knowledge on either ||Va;||£2 ^(^jj ) or \\Va\\i,2 ^(Hu)- 
But they are controlled in next section when we derive estimates on V and V^a. 
Proof. We only sketch one estimate (which may be regarded as the most difficult one). 



The rest is similar so we skip. For Vsa, according to (2.10), we have 



Vga = -^trx« + V®/3 + Aua - 3(xp +*X(^) + (C + 4r/)§/3 
= trx, ■ a + V/3 + ■ a + V ■ ^0, 



The term ij^g ■ denotes a sum of products of this form. We recall once again that 1^ denote a 
generic connection coefficient and ^ denote a generic curvature component. The subindex g means this 
component is not anomalous. 
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We can bound V/3 by TZ^ and bound 



For trx^a, in view of Proposition 2.5 
above, we have the desired estimates 



ast two terms by C6i thanks to Holder's inequahty. 
, it is bounded by Tq6^2 + C6i. Combining all the 
for Vsa. This completes the proof. □ 



Combining previous two lemmas, we derive the comparison estimates for null compo- 
nents of DsRa/B-yS and D^Ra/s-ys- 

Proposition 3.9. (1). For a{D^R), we have 

[sc) ^ {sc) ^ — ^i'' 

(2). For^iDj^R) ^ {a{D^R),a{DiR),a{DsR)} , except l3{D^R) on and (3{D^R) on 
Hu, we have 

||vi/(D^i?)IU.^^^(^^) <7^o + 7^l + C(5^ \\^>{DNR)\\L2^^^^nj<n, + n, + C5-^. 

We now consider horizontal derivatives and compare Vc^l/ with ^[DcR)- 
Proposition 3.10. If^^(3, we have 

\\^{d,r)-vMlj r^,) + ||^pci?)-Ve^lk2 ,^^<n + n + C6K 

mD,R)\\^. (^) + ||vl>(D,i?)|L. <7^ + ^ + C5i. 



Proof. The proof is quite similar to Lemma 3.6 Notice that, if \l/ 7^ /3, then 

^{D,R) - Ve^ = ■ ^ + trx^, ■ ^g. 

The product term is bounded by C5^ as before and last term is bounded by 7^ + ^ by 
definition. Thus we have the first inequality. The second follows immediately. □ 

We turn to the anomalous term Vc/3. First of all, notice that 

(3{DcR)a = ^cl^a - ]^X^iP^ba + Cc ' 1^ a - ]pCca ' (3p - O") 

= trxg ■ a + V/3 + ^3 ■ ^ + tr^^ ■ 
Thus, the first term trx^ ■ a causes anomalous behavior. Hence we have 

W{d,r)\\l^ <h) + W{d,r)\\l. , <i,.5~'^+n + n + csl (a.i) 

3.3. Mild Anomalies. We see that a{D4R), a{D^R), a{D^R) and (3{DcR) are anomalies 
since their estimates cause a loss of 6~2. We shall distinguish a^D^R) and f3{DcR) from 
a{D4^R) and a^D-^R). Heuristically, a{D^R) and f3{DcR) are anomalous in zero order; 
alD^R) and a{D^R) are anomalous in first order. 

We observe that alD^R) and P{DcR) are anomalous because they can be written 
schematically as 

trx, ■a + ^g + (V^), + ■ ^. (3.2) 
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This amounts to say that their anomahes come from the zero order term trx^ ■ a. 



We 



shall call a{DsR) and l3{DcR) mild anomalies. By the above comparison lemmas, Vaa is 
also a mild anomaly. 

In what follows, we use this observation to replace a^D^R) and (3{DcR) by the worst 
term a. We have two advantages of performing this substitution: the first is that we have 
already controlled a in Proposition 2.5 the second is that we can use L^^^.^ estimates on a 



(but not on either a{DsR) or P{DcR)j. The other terms are much easier to control in all 
the contexts since they are not anomalous at all. This becomes much clear and important 
when we derive estimates in following sections. 

4. Anomalous Estimates 
Y 



4.1. Estimates on V^a. We take X 
derive 



Z = L in (2.25) and = 64 in (2.26) to 



^ \a{D,R)\^+ [ \I3{D^R)\^<[ \a{D^R)\^ + [[ D^Q[D^RU^^ + '^'^\^,Q^\^. 

Hu Jh^ JHo JJv(u,u) 



In view of Lemma 3.6, by ignoring the term ^{D^R) on the left hand side and passing to 



scale invariant norms, we derive 



l^Mll, .„(o,.) < ||V4a||' ,^(0,.), + C5^+5' ^ D'^g[D4i?]M44 + ^%.Q^%4 

(sc)*- " (3c)( ) JJt>(u,u) 



' JJV(u,u) 

< b-^Zl + {I + J) + K. 

We split the bulk integral into three terms J, J and K in such a way that the integrands 
of / + J are exactly the terms appearing in D^Q[D4/2]^444 and the integrands of K are 



exactly the terms appearing in ^^^ti^uQ 



44. 



Remark 4.1. In the above terms I, J and K , there is a 6 to certain power in front of 
the integral which reflects the fact that the inequality is scale invariant. In what follows, 
this phenomenon appears constantly. To simplify expressions, we shall always ignore this 
6 to certain power since they can be easily restored by investigating the signatures. 



According to (2.23) and ( 2.26[ ), the integrand in I + J has the following schematic form, 
L'4i?4^"4tt + D^m/i" *J^^l = ^{D^RY'"^ ■ ■ ¥"'^ + ■ <i!{D^RY''^ ■ ^{D^RY''^ 

with So + si + S2 = 6. Here we use the upper index Sj's to indicate the signature. We 
require the integrand of I is of the form \l/(D4i?)'-^°) ■ vl>(*i) . \l/(*2) g^j^^ ^j^g integrand of J 
is of the form ■ "^{DiR^'^^ ■ '^{D^R)^'^\ We remark that ^'■'"^ ^ trx^. F 

We now examine the possible anomalies in the integrand of I. For '^(D^R)^^'^^ on H^, 
according to Lemma 3.9, the only possible anomaly is alD^R) which has 3; for "^^^ and 
\E'**^, the only possible anomaly is a with signature 2. Thus, in view of the signature 



^■^ In following sections, we shall always split the error integrals into three terms /, J and K in the 
same way. And we shall not explain when one encounters this situation later. 
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consideration, at least one of the terms is not anomalous. We say that, in this situation, 
we do not have triple anom alies. 

In view of Corollary 3.4, we estimate / as 

followf3 

J 

Jo Jo 
For K, according to signature considerations, we can also avoid triple anomalies. Thus, 

Q^\4 = ^ ■ aiD^R) ■ a{DiR) + ■ ^g^D^R) ■ ^{D^R). 

We dominate ipg and ip in L^^^ norm to derive 

Jo Jo 
For J, it possesses most of the features of K. We still need to investigate the structure 
of the integrand. In fact, the integrand is of the following formp^ 

D^R/4,^ ■ D'^L'' ■ D^R^p^s = ^^'"^ ■ ^(1^4^?)^''^ ■ "^{D^RY'^^ (4.1) 

Thus, since ip^^o^ ^ ^^^o' bound J by 

j<C6'^ .H„)<C55 / criers <cri 

Jo Jo 
Combining the estimates for J, J, K and Young's inequality, we derive 

l|V4«f^. ,^,0,.,. < r%^ + cri < 6-'c{io) + C6l 

In conclusion, we have the desired anomalous estimates on W^a 

.Mo,^)< 5-h{Io) + C5-K (4.2) 

(SC) V ^ / 

Remark 4.2. T/ie a6ofe analysis also yields the following estimates which is absent in 



Proposition 3.9 



^(0,.) <r5Xo + cr. (4.3) 



^'^ In this section, since we deal with exceptional terms, we can be extremely wasteful. In other 
sections, to derive energy estimates, we have to exploit the cancelation hidden in the integrand. 

We ignore the terms coming from the Hodge dual of curvatures since due to the symmetry. They 
can be treated exactly in the same way. This remark applies to all situations with one exception at the 
end of the paper. 
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4.2. Estimates on Vsa. We take iV = L in ( |2.26[ ) and X = Y = Z = L in (\2.25l to 
derive 



fiu^ 33- 



/ \PiD;R)\'+ [ \a{D,R)\^<f \^iD,R)\^+ [[ D'^Q[DsRUs^ + ^^\,,Q 

JHu JH_u J Hq JJv{u,u) 

In view of Lemma 3.6, by ignoring the term ^{D^R) on the left hand side and passing to 
scale invariant norms, we derive 

l|V3«||' .^(0,.). < r% + C6'^ + 6-' [[ D'^QlDsRU.^ + ^-VQ^"33 
<6-^Xl + I + J + K. 

We split the bulk integral into J, J and K in the obvious way as before. The integrands 
in / have the following schematic form, 

with So + Si + S2 = 1. According to signature considerations, none of ^(^i) and ^(^^3) can 
be anomalous since the signature for anomalous curvature component has signature 2. 
Thus, we estimate I as follows 

I<6h-'£mDsRY-nLl^^^n^^^^^^^^ 

ru 

<sh-^ / CS-'^C^du' <C. 



We now deal with K. Again, by counting signatures, the anomalous terms ai^D^R), 
ai^D^R) and f3{DcR) do not appear in the integrand since their signatures are at least 
2. Thus the only possible anomaly is a{D^R). Hence, if their is a double anomaly in 
curvature components, they must be two a{D^R). This forces the connection coefficient 
to have signature 1. In particular, it can not be trx^. Thus, the integrands have the 
following schematic form 

Hence, bounding ip in L^^^ yields. 



K < (CTi + 1)5-' rii^'(D3i?)iu^^^,(i.jii^'.(i^3i?)ik^^^,(,.j 

t/ 

+ C5'^ ■ r ||a(D3i?)IU^^^,(^j||a(^3i?)IU^^^,(Hj < Cr ^ 



The analysis for J is exactly the same as K, hence, we have 



Recall that trv is a constant with worse estimates in Lf . than usual connection coefficients. 
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Combining those estimates for /, J and as we have done for V4Q; in last subsection, 
we derive 

||V3a||^2 .^(0,.) <r^Xo + C53. (4.4) 
Remark 4.3. As a byproduct, we also have the following estimates which is absent in 



Proposition 3.9 



||V3^||^2 .^(0,,) <r5Xo + CT3. (4.5) 

5. Energy Estimates on Outgoing Derivatives 

In this section, we derive energy estimates for V4\l/. We first take N = L and all 
possible X,Y,Z E except for the case {X,Y,Z) = {L,L,L), in (2.25) to deriv^ 

|^(Ai?)^'^|'+ / |^(Ai?)^'"5)|2 = f |^p^/^)W|2 



+ 



JJV(um) 



We then multiply both sides by 6"^^^^ and sum all those inequalities, in view of Remark 



4.1[ we have 

J2[\\^{D,R)\\l, .^(o..K + ||^(i^4i?)||J. (o,.K]=X2 + / + J + ir (5.1) 

(sc) ^ ' (sc) ^ — ^ ^ 

where 



I + J= [[ D>^Q[D,RlxYZ, K= [[ ^^'\,,Q[D,R] 

JJv{u,u) JJv{u,u) 



IIU 

NiN2- 



The sum in (5.1 ) is over on all possible ^iD^^RV?, except aiD^R) and W/SiD^R) 11 2 /rr(o.")v 
The definition for J, J is in an obvious way as in the previous section. 

5.1. Estimates for /. The integrand of / can be written schematically as 

with total signature 3<so + si + S2<5. We shall further split / into 

I = I0 + h + I2 + h 

where Ik denotes the collection of the terms whose integrand has exactly k anomalies. 

We observe that J3 = according to signature considerations. Otherwise, for a triple 
anomaly, thus Sq > 2, si = 2 and S2 = 2 which contradicts the fact that Sq + Si + S2 < 5. 
We also notice that the only possibility for "^(D^R)^'^"^ to be anomalous is "^(D^R)^^'^^ = 
a{DiR). 
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We shall use letters Nq, Ni and N2 to denote either L or L coming from X, Y and Z. 
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5.1.1. Estimates for Iq. The control of Jq is much easier than other terms. 

lo < 5^ £ IMD^Rt^^h^^^^^^^^^^^^ < C5l (5.2) 

5.1.2. Estimates for Ii. We further spht Ji into two terms 

Jl = 111 + /i2 

where the corresponding integrands for In and lu are 

^(D4i?)J'°^ ■ a ■ ^^''^ and a{DiR) ■ ■ ¥g'^^ 

respectively. 

111 can be bounded as follows, 

112 requires an extra integration by parts. Before going into details, we first show how 
to replace a{DiR) by V ^a. According to the proof of Lemma 3.6 , a{DiR) — V 4,a = ipg-'^ . 
Thus, 

/l2 = // V4a ■ ■ ^[''^ + [[ ■ ^ ■ ^^'^^ ■ = h21 + /l22 

JJV(u,u) JJv{u,u) 

For Ji22, since the nonlinearity gains one more term, we gain an extra 6^ thanks to scale 
invariant Holder's inequality. Thus we can bound it as follows. 



2 



Remark 5.1. The above analysis allows us to freely switch between \l/(Djv-R) one? Vat^I/. 
The difference contributes at most C5^ to energy estimates. In the rest of the paper, we 
shall ignore this C6^ and we shall not distinguish "^IDnR) from Vn"^- 

We move on to lu. 

/i2 = /i2i = // V,a-^g-^g= II a- V,^g ' g + l " " " ^9 l^^o 

JJV(u,u) JJV{u,u) JH^, 



For the bulk integral, in view of Lemma 3.7, we bound it as follows 

a ■ V,^g ■ < 5i / ||V4VI/,|L2 (^^,)||a|L4 (^^,)||vl/,|L4 (^^,) < CS'^. 

'©(Ujij) Jo 

For boundary integrals, we bound them as follows 

Thus, we conclude I12 ^ C5^. Combining the estimates for In and /12, we derive 

/i<C5i. (5.3) 
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5.1.3. Estimates for 12- We split I2 into two terms 

I2 = I21 + -^22 

where the corresponding integrands for I21 and I22 axe 

*(L>4i?)J^o) -a- a, aiD^R) ■ a ■ ^^'^^ 

respectively. 
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We first show that hi = 0. In view of ( |2.23| ), ( |2.26[ ), ( |2.27[ ) and ( |2.28[ ), the integrand 
of /21 has the following expression 



DnRx^z^ ■ R^ny'' ■ R, 



(5.4) 



We write only one term in (5.4). But this term is generic and the rest can be treated in 



the same way. If we look for a ■ a, none of or z/ can be 3 or 4; this forces Y = (3 = 6 = 4. 



Thus (5.4) reduces to 



DaRxszs ■ a ■ a. 

If X = L or y = L, it vanishes; if X = Y = L, the the total signature is 7, it is impossible. 
Hence I21 must vanish. 

For I22, since total signature is at most 5, thus S2 = 0. Thus the integrand of I22 niust 
be EH 

a{D4^R) ■ a ■ a. 



In view of Remark |5.1[ we have 

1 



'22 



V4Q; ■ a ■ a 



V{u,u) 



VJa ■ a) ■ a 



V{u,u) 



a ■ a ■ V^a + a ■ a ■ a 



V{u,u) 



\u'=0 



For the bulk integral, we use Bianchi equation (|2.19|) to replace V^a. One can argue as 
in Remark 



5.1 



only the top order term V®/3 in (2.19) is relevant. The quadratic terms 



in (2.19) contribute at most a C^*. By performing an extra integration by parts, we can 



move V to a to eliminate the anomaly of a. More precisely, we proceed as follows 

a-Va- 13\+ C53 



a- a- V4a| < I // " ■ " ■ VO/3| + 05^ < 

V{u,u) JJV(u,u) 



V{u,u) 



< 6^ 



iVal 



To control the boundary integral, we need to improve Proposition 3.2 
Lemma 5.2. If 6 is sufficiently small, then 

\MLf AS) < 5-3c(Xo) + r 3c(Xo)7^^ + 



17 



Since a and a are traceless, ■ a ■ a — ^ ^cia ctb'^ a J^. 
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Proof. In view of Proposition 3.3, we have 



Combined with Proposition 2.5 and (4.2), thus 

< (^0 + C6-^ + ||v«|U2^^^(j,„))5((i + r 5)Xo + c<53)5. (5.5) 

We then bound ||Va||i2 (j:^^) by 7^ to complete the proof. □ 

(sc) 



Thus, combined with Proposition |2.5[ the boundary integral is controlled as follows 

< 55 (Jo + c(Xo)7^^ + 3c(Jo) + r 3c(Xo)7e5 + C6^^y 

<c(Xo)(l + 7^i) + C5i. 
Thus, we have the following estimates for I22 and I2, 

/2 = /22<c(Xo)(l + 7^t)+C^i 

5.1.4. Conclusion. Combining the estimates for Ji and I2, we derive 

J<c(Xo)(l + 7^i) + CTi. (5.6) 
5.2. Estimates for K. The integrand in K can be written schematically as 

with total signature 3 < sq + si + S2 < 5. The appearance of trx^ is due to ^-^tt. We split 
K into 

K = K0 + K1 + K2 

where K^. denotes the collection of terms with exactly k a{D4R)^s. We observe that 
K2 = according to the signature considerations. 

5.2.1. Estimates for Kq. We first split Kq into two terms 

K2 = Koi + K02 

where the corresponding integrands for Kqi and Kq2 are 

^(^0) . . and trxg ■ '^{D^R)^''^ ■ "^{D^RYg^^ 

respectively. 

Kqi can be bounded as follows. 
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Kq2 can be bounded in a similar way, although we lose a due to trx^, 

J 

r\\^{D,R), 

Jo 



Once we plug Kq2 into ( |5.1[ ), in view of Gronwall's inequality, it can be absorbed by the 
left hand side of (5.1). Hence we can ignore K02 at this stage and simply write 

Ko<CS"^. (5.7) 

5.2.2. Estimates for Ki. First of all, we remark that trx^ does not appear as connection 
coefficients in the integrand of Ki. Otherwise, in view of (2.24), the appearance of trx^ 
is through 

Since a[D^R) appears, this forces Ni = N2 = L. Thus, the integrand can be written as 

■ QiD.Rr",, = x^\paab - 2a*aab). 
Since a and *aab are traceless, the trx^ will be canceled. This argument also shows 

X ■ a{D4^R) ■ '^{D^R^g'^^ is absent in the integrand, otherwise they must be either x " 
a{D4^R) ■ pi^D^^R) or x ' cn^D^R) ■ a{D4^R) with total signature 6. This is impossible. 
We further split Ki into two terms 

= Kn + Ku 

where the corresponding integrands for Ku and K12 are 

X ■ a{DiR) ■ ^{D^RYg^^ and ^(""^ ■ aiD^R) ■ ^{D^RYg^^ 

respectively. 

For Kn, as we just observed, '^{DiR)'^^'^ = p^D^R) or ^{D4R)'g ''^ = a^D^R). We only 
treat the first possibility. The second one can be estimated exactly in the same way. In 
view of the Remark 5.1 and (2.15), we have 

Kii= // x-V4a-V4P= // X- V4a(-^trxp + div/3- ^x-a + C-/3 + 2r7-/3) 

JJv{u,u) JJv{u,u) ^ ^ 



X-V4a-V/3+ // {x-^'ia)ix-a) + C6-^ = Kui + Kii2 + C6-K 

V{u,u) JJv{u,u) 

For Kill, we reduce the number of anomalies by an extra integration by parts. 



/iii=// V4X-«-V/3+// x-«-V4V/3+/ x-a-V/3|5l^ 

'V{u,u) JJv(u,u) Jh^, 



Klin + i^iii2 + K 



1113- 
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To estimate Kuu, we use null structure equation (2.8), 
Kuu= II ( V% + 2a;x - Jtrxx + !Z% - ^t^^XX) ■ a ■ V/3 = ff x ■ « ■ V/3 + 

JJv{u,u) ^ ^ jJv{u,u) 



IV{u,u) JJv{u,u) 

To estimate Kuu, we compute the commutator 



We claim, 



[V4, V]/3 = -X ■ V/3 + */3 ■ /3 + ^ (r/ + r/) V4/3 + ■ /3 ■ X, 
I // X-a-[V4,V]/3|<C53. 



In fact, the worst term should be *(3 ■ /3. We can use L?° ^ norm on x, Lf . norm on two 



(sc) ^^^^^^^ A' ^{sc) 

/3's and L^^^^ on a to dominate it. Thus, up to an error of size CS^, 



Kui2= // x-a-VV4/3= // Vx-a-V4/3 + // x ■ Va ■ V4/3. 

JJv{u,u) JJv{u,u) JJV{u,u) 

For the first integral, 

jjt>{u,u) Jo 
For the second integral, 



JIvi 



X ■ Va ■ V4/3 < 5^ /"l|V4/3|U^^^^(H„,)l|Va|U.^^^(H^,)||x||L^^,(//„,) < 



Hence, 

i^lll2 < C5 

To estimate -ft'ina, we proceed as follows. 



i^iiis < 



X ■ a ■ V/3| < 5^||V/3iU2^^^(^j||a|L4^^^(^j||xlL*^^,(|/j 



< 6-^n{6—^c{io) + r3c(Xo)7^2 + c5T5)(r j + c^?) < c(Xo)(7^ + n)^ + cs-i. 

Combining the estimates for Kim, Kuu and Km^, we derive 

Kui<c{Io)nl+C6l 

We now turn to Kiu- 
Ku2= II (X- V4a)(x-«) = - ff^ (V4X-a)(x-«) + J // V4|x-«r 

JJv{u,u) JJv{u,u) JJV{u,u) 

(V4X ■ • «) + ^ / V4|x-«r ly'Io = ^1121 +^1122- 
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To estimate -ft'1121, we use (2.8) to derive 
^1121 = [[ (V§r7 + 2ujx - Jtrxx + rpV - ^t^xx) ■ a ■ X' a = [[ + C6^ 

JJv{u,u) ^ ^ JJV(u,u) 

To estimate -ft'1122, we proceed as follows 

Thus, we have K\\2 ^ Cb^. Hence, 

irll<c(Xo)(7^ + R)i+C5l 

It remains to control K\2 whose integrand is tll^g^^ ■ alD^R) ■ "^ID^R^I^^K The key[^ is 
to observe ipg 7^ u. Otherwise, since total signature is at most 5, it forces the integrand 
to be a; ■ a{D4^R) ■ a^D^R). Once again, we investigate the origin of this term 

By direction computation, we see that ai^D^^R) -a^D^R) never appears in Qi^D^RY^ m^n^, 
hence we arrive at a contradiction. 

Based on the above observation, since we have L'^^sc) estimates on V^ipg, we can proceed 
exactly as we did for Ku to derive 

In this case, since ipg is not anomalous, the estimates are much easier to derive. We skip 
details. Combining the estimates for Ku and ^^12, we derive 

Kl<c(Xo)(7^ + R)i+C(53. (5.8) 
5.2.3. Conclusion. Combining the estimates for Kq and i^'i, we derive 

if < c(Xo)(7^ + E:)^ (5.9) 
5.3. Estimates for J. The integrand of J can be written schematically as 

with total signature 3 < sq + si + S2 < 5. We split J into 

J = Jo + Ji + J2 

where Jk denotes the collection of terms with exactly k anomalous curvature components. 

We make the following observation: the only anomaly for ^![DiR)'^^'^^ is a{DiR)] things 
are different for ^!{Di,R)^^'^\ it has multiple choices: a{D4R), a{D^R), a{D2,R) and 
P{DcR). We also recall Section 3.3 which asserts that last two anomalies are mild in 

In [U, the authors derived all -^^g^) estimates for V4V' except for the case ip — tu. 
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the sense that they all come from zero order term a. This will be extremely crucial in 
our proof. 



5.3.1. Estimates for Jq. The integrand of Jq, by definition, has the following schematic 
expressions 



Thus, we bound Jo as follows, 

J 

5.3.2. Estimates for Ji. We split Ji into Ji = Jn + J12. For Jn, ^{D^R)^'^^ = a{DiR); 
for J12, ^HiD^RY'^^ = '^{DiR)^"'^ is not anomalous. 

The integrand of Jn has the following form (recall that we ignore the Hodge dual part) 

DiRx'^Y^ ■ D^V ■ DyR^pzs. 

Since ^{D^R)^'^^ = aiD^R), it forces X = Y = L, Z = L, p = a and 6 = b. Thus, the 
integrand is reduced to 

V4a ■ D^L" ■ D,R^a3b. 

Since D'^L^ = 0, we see u ^ 3. If = 4, we use Bianchi identities, up to lower order terms, 
to convert D^R^aSb into the form V^E'g. Hence, the only possible forms of the integrand 
are either x " ■ V\E'g or ipg ■ V4,a ■ V\E'c, {^ipg 7^ u). We then split Jn as 



Ju = x-V^a- + // ^g- Via ■ V^'^ = Jm + J112. 

JJv{u,u) JJv{u,u) 

Notice that J112 has exactly the same form as K12 term in the previous section, thus 



The estimates for Jm is quite similar to Kn. First of all, according to signature consid- 
eration, V\E'g = V/3, thus 

Jin =11 V,x-a-V^+ [[ x-a- [ x-0i-Vf3_ 

JJv{u,u) JJv{u,u) J H_^i 

= Jiiii + J1112 + JlllS- 



To estimate Jm, we use null structure equation (2.1), 

Jiiii= If ViX-a-V§_= jj (-trxx-2a;x-«) 

JJv(u,u) JJv(u,u) 

a-a-V/3 + CS^ = // a-Va- (3 + C5^. 



'^{u^u) JJ'V(u,u) 

Hence, 
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To estimate J1112, we compute commutator 

[V4,V]^=-X-V^ + *^-^+^(?7 + 77)V4^ + 77-^-X, 



Using this formula, we can easily derive 

I // x-«-[V4,V]^| <C5i 

JJv(u,u) 

Hence, up to an error of size CS^, we have 



Jiii2= // x-a-VV4^= // Vx-a-V4^+ // x • Va ■ V4^ < C^J. 

JjV{u,u) jJV{u,u) JJv{u,u) 

To estimate Jma, we proceed as follows, 

Jni3 < I / X • « • V^l < 5^||V^|U2^^^(^j||a|U4^^^(^j||x||L4^^^(^^j 

< 55^(r 3c(Xo) + r 3c(Xo)7^^ + C5^)(r 3 + C53) < c(Xo)(7^ + + c^j. 

Putting things together, we have 

Jn<c(io)(n + n)-^+c5l 

We move on to the estimates for J12 where ^{D,yR)^^^^ must be anomalous. Thus, 
because D^'L^ = 0, it must be either a{D4^R) or (5{DaR). We can assume ^ 
a(D4R). Otherwise, by signature considerations, those terms (of the form V\&g ■ ■ V4Q;) 
have already been treated in Jn, so we ignore them at this stage. The anomalous term 
must be (3{DaR). As we emphasized at the beginning of the section, its anomaly comes 
from the lower derivative term a. Hence, up to an error of size (7^4, we can replace the 
integrand by 

where V^g 7^ Va according to signature considerations. The second term can be esti- 
mated easily by placing L^^^^ estimates on ipg and a (which saves an extra 5*), thus we 
bound J12 as follows 

J12 = C53 + // ip-a- V^g = + // ■ a ■ + // ^ • Va • *g 

JJv{u,u) JJ V{u,u) JJv(u,u) 

«/ 

Finally, we can bound Ji by 

^i<(^ + ^)^+C5j. 
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5.3.3. Vanishing for J2. We claim J2 = 0. Since both curvature terms in the integrand 
are anomalous, it must be 

D^R^'J' ■ D^L" ■ D^R^aSb, 
Since total signature is 5 at this situation, the second anomaly must be either alD^R) 
or P{DcR). Their signature must be 2. Once again, since D^Li^ = 0, thus the second 
anomaly must be P{DcR). Thus, u = c. Finally, since the signature of DcRfj,a3b is at most 
|, it can not be anomalous. This shows J2 = 0. 

5.3.4. Conclusion. Combining the estimates for Jq and Ji, we finally derive estimates for 
J as follows, 

J< c(Xo)(7^ + ^)i +C5i. (5.10) 
We combine (5.1), (5.6), (5.9) and (5.10), we conclude that 

J2miD,R)\\l, ,^,,.u + \\^iD,R)\\l, ]<c(Xo)(l + (7^ + E)^) + C^5i (5.11) 



Considering ||/3(Z)4i?)||^2 (/i-(o>")) (5.11), we have 



{sc)>.— It 

|V4^|L, ,^(0..) <c(Xo)(l + (7^ + ^)i) + C(55. 



In view of (2. 17), modulo quadratic terms which gives an error of size C6i, we have 
II - Vp +*V(t| 



< 



c(Xo)(l + (7^ + E:)5) + C5 



According to the standard elliptic estimates for Hodge systems, we derive 
II (Vp 

Considering \\l3{DiR)\\ 



||(Vp, Va)||^, (0,.) < c(Xo)(l + (7^ + ^)5) + C5s. 



n). in (5.11), similarly, we derive 



||(Vp, Va)||^, (^(0,.)) < c(Xo)(l + (7^ + ^)5) + C6-S. 



Considering ||/3(D4i?)|| 2 (5.11) and (2.11), we derive 

II Va|| 2 ,^(0,.). < c(Xo)(l + (7^ + E:)5) + 

(sc.) V / 



Considering ||p(Z}4i?) II 2 frrC),")-,, ||c"(-D4i?) IL2 ,„(o,«)-, in (|5.11|) and (|2.15|) as well as (|2.13|), 
we derive 

Similarly, we have 



l|V/3||^. .^(0,.). < c(Xo)(l + (7^ + R)i) + C6"s. 



<c(Xo)(l + (7^ + ^)2) + C5B. 
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Putting all them all together, we have 

||(Va,V/3,Vp,Va)|L. + ||(V/3, Vp,Va)|L. 

<c(Io)(l + (7^ + R)i) + C(5i. (5.12) 
6. Energy Estimates on Incoming Derivatives 



We take N = L, (X, Y, Z) = (L, L, L) or (X, F, Z) = (L, L, L), in ([2^ to derive 



JJvU.u] 



iV(u,u) 



We multiply both sides by 6'^'^ ^ and sum those two inequalities, in view of Remark 
we have 

mD,R),piD,R),aiD,R))r + IMD.R), aiD.R), PiD,R))f (o,.). (6.1) 

<Xl + I + J + K 
where 

/ + J = I // D^'Q[D:,R]^xYzl K= \ II W'V^,Q[D3i?]^''jv,7v,|. 

JJV(u,u) JJV{u,u) 

In this section, the possible anomalies for ^!{D^RY^^'> are a{DiR) and a{D^R)] the 
possible anomalies for ^(Dj,i?)(^2) ^re a{DiR), a{D-iR), a{DiR) and (3{DcR). We still 
emphasize that the anomalies for a{D^R) and P{DcR) are mild in nature. They come 
from the lower derivative anomaly a. For example, in applications, we always use the 
following expressions, 

aiD^R) = a + C6K PiD^R) = a + C^^ . (6.2) 

6.1. Estimates for /. The integrand of / can be written schematically as 

'^{DsR)^''^^ ■ "^^'^^ ■ \E'^''2^ 
with total signature sq + si + S2 = 2 or 3. We split / into 

J = Jo + /i + /2 + /a 

where Ik denotes the collection of terms with exactly k anomalies. According to the 
signature considerations, J3 = 0, otherwise S2 + S3 > 4. 

6.1.1. Estimates for Iq. We control Iq as follows, 

h < 5^r^ j mD^Rt^Ll^^^H^^^^^^^^ < C5l. (6.3) 

J 
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6.1.2. Estimates for Ii. We split Ji into three terms 

h = hi + h2 + -^13 
where the corresponding integrands for In, I12 and /13 are 

respectively. 

To estimate In, we have 



To estimate J12, in view of (6.2), modulo C5i, we have 



To estimate J13, we observe that none of the \E'g's is a. Otherwise, since the total 
signature is at least 2, this force another g to be a which is impossible(in this case, we 
only assume Ji contains one anomaly). Thus, up to an error of size C5^, 



h,= V,a-^g-^g= a-Vs^g-^g+ a-^,-^,\:,To. 

JJV{u,u) JJv{u,u) J H^i 

For the bulk integral, we bound it as follows 
For boundary integrals, we bound them as follows 

1 



Hence, we bound /13 < C5^. Combined with the estimates for Jn and I12, we derive 

/i<C5i. (6.4) 



6.1.3. Estimates for h- According to signature considerations, the integrand of I2 must 
have the following schematic form 

aiDsR) ■ a ■ 

Since the total signature is 2 or 3, thus \E'c, G {a,p,a}. The estimates for = p and 
\E'g = (T are similar, we shall only concentrate on \E'g = p. Hence, we further split I2 into 



h = hi + hi = I Vsa ■ a ■ a+ // Vsa ■ a ■ p. 

'v(u,u) JJV{u,u) 
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To control I21, we use move the derivative to a better position, 

hi = // V3(a ■«)■«= // a-a-V3a+ / a-a-a|"'r" 



V{u,u) JJV{u,u) J H^i 



lii'=0 



In view of (6.2), up to an error of size C5^, the bulk integral is 



a-a-a 

. 1 



Then we can place Lf norm on a to save an 5-1 to derive 



hi^A II a-a-a\ + \ I a-a-a\<C6*. 



To control I22, we still integrate by parts. 



h2= // V3a-a-p= // a-p-V3a+ // a ■ Vap ■ a + / a-p-a|Il'=o- 

JJv{u,u) JJV{u,u) JJV{u,u) J H^i 

Exactly as for /21, we can bound each term by C5^ easily. Thus, 

h<hi + h2<C5-^. (6.5) 

6.1.4. Conclusion. Combining the estimates for Ji and J2, we derive 

I<C5-^. (6.6) 

6.2. Estimates for K. The integrand of K can be written schematically as 

(^ + trxg)^'"^ • ^(/^3^)^''^ ■ ^(1^3^)^''^ 
with total signature sq + si + S2 = 2 or 3. We split K into 

K = Ko + Ki + K2 

where Kk denotes the collection of terms with exactly k anomalous curvature components. 

6.2.1. Estimates for Kq. We split Ko into two terms 

K2 = Koi + K02 
where the corresponding integrands for Koi and Kq2 are 

■ vl/(D3i?)(^i) . ^{DsRY;'\ trx„ ■ ^{DsRY;^^ ■ vl/(D3i?)^^ 

respectively. 

Kqi can be bounded as follows. 



Koi < 52 
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Ko2 can be bounded as follows, 



(S2)|2 

a I 



Ko2 < I // trv ■ vl/(D3i?)(«i) . ^{DsRY;^^] < [[ |vl>(D3i?)^f + I^Psi?) 

JJv(u.u) JJv{u.u) 
pu pu 

According to Gronwall's inequality, Kq2 can be absorbed by the left hand side of the sum 
of (5.1) and (6.1). Thus, we can ignore /02 at this stage. Hence, 

Ko<C5K (6.7) 

trx, 



it Ki into 



6.2.2. Estimates for Ki. One can argue exactly as in the beginning of Section 
does not appear as the connection coefficients in the integrand of Ki. We sp 
two terms 

K^ = Ku + K,2 
where the corresponding integrands for Kn and K12 are 

respectively. 

We first control Ku. Up to an error of size CS*, using Bianchi equations, we have 
Kn= [[ ^ ■ « ■ (Vs^)^^^) = [[ ^-a-V^g 

JJv{u,u) JJv{u,u) 



V{u,u) JJv{u,u) 

Now it is routine to bound those two integrals by C5^. Hence, 
We now control K12, we need to further decompose it into 

K12 = Kui + K,22 =11 X- cy{DsR) ■ ^{DsR)g + [[ ^- a{D^R) ■ ^{D^R)g 

JJv{u,u) JJv(u,u) 

where ip ^ X- We remark that x does not appear, otherwise ^{D^R)g = alD^R) which 
leads to a double anomaly. 

To estimate J121, observe that the appearance of x is through the term 



NiN2- 



Since a^D^R) also appears, this forces Ni = N2 = L. Direct computations show 
^{D^R)g = p[D^R) or '^{D^R)g = a{D^R). Without loss of generality, we may assume 



K121 = I X- Vsa ■ Vsp. 
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We use (2.16) to replace Vsp. Up to an error of size C^^, we have 



^121 = // X ■ ■ V/3 = /i211 + /l212 + ^1213 

'V{u,u) 



V{u,u) JJv{u,u) Jh^i 



u =u 



For the boundary integral /1213, 



^1213 < sup| / x-a- V/3| < 5^||xlU-)l|tt||L2^^j(/f„)l|V/3|L2^^^(^^^) < C^i 
For -ft'1212, we compute commutator 

[V3,V]^=-X- V^ + *^-^+^(^ + r2)V3^ + X-^-^- 

All terms are quadratic except the first one which is essentially linear in V/3. Thus, it is 
routine to derive 

I [[ X-a-[V3,V]^|<| [[ x-«-V^|+C52<C5i. 

JJv{u,u) jJv(u,u) 

Hence, up to an error of size C^J, we have 



^1212=// X-a-VV3/3= // Vx-a-V3^+ // X ■ Va ■ V3^. 

JJv{u,u) JJv{u,u) JJv{u,u) 

Hence, we eliminate all the anomalies. Thus, it is routine to bound both integrals by C5^ 
as well as 



For -ft'1211, we use (2.9) to replace V4X, thus, 

^1211 = 11 (-^trxx + V®r7 + 2ux - Jtr^X + V^V-) ■ a ■ V/3 

JJv{u,u) ^ ^ 
lV{u,u) JJv{u,u) 



2 -0 



For the first integral, we simply use L^^-j on x fo overcome the anomaly; for the second 



one, we control it as follows. 



jjt>{u,u) Jo 

This yields the control of lun- Combined with the estimates for -ft'1212 and -ft'1213, we 
have 

To estimate -ft'122, we observe that in its integrand ip ■ a{D^R) - ^ {D^R) we have ip^to. 
Otherwise the total signature is less than 2. This allow Vs'?/' to have L^^^^ estimates (while 
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V3U does not, see [3]). We can proceed exactly (much easier in reality) as we did for K121 
to show that K12 < C6i. Combined with the estimates for Kn, we derive 

6.2.3. Vanishing for K2- We claim K2 = 0. The integrand of K2 may have the following 
three types of terms 

{ij+tTX^)-a{D,R)-a{D3R), {ij+iix^)-a{D-,R)-a{D,R), i^+tix^) ■a{D,R)-a{D,R). 

Since the total signature is either 2 or 3, this rules out the first two cases immediately. 
To rule out the third one, notice that its appearance is originally through the following 
term 

^'''^TT^.QlDsRrN.N, 

Here, A^^i and N2 are null. Once again, direct computations show that a ■ a never appears 

in Q^uNiN2- 

Finally, we conclude that 

K<C5-K (6.8) 
6.3. Estimates for J. The integrand of J can be written schematically as 

(V- + trXg)^'"^ ■ ^{D^R)^'^'^ ■ ^{D^Rf'^'^ 
with total signature sq + Si + S2 = 2 or 3. We split J into 

t/ = t/o + t/i + t/2 

where Jk denotes the collection of terms with exactly k anomalous curvature components. 

6.3.1. Estimates for Jq. The integrand of Jo must be one of the following forms 

v-^^") • ^(/^3i?)^^ ■ ^{D,RYp\ trxg ■ ^(^3^?)^^ ■ "^{d^rY;'^. 

Thus, those terms can be dealt with exactly as we have done for Kq is last subsection 
(of the current section). Thus, modulo the terms which will be removed eventually by 
Gronwall's inequality, we derive 

Jo<C6l (6.9) 

6.3.2. Estimates for Ji. According to the position of the anomaly, we split Ji into 

Ji = Jii + J12 + J13, 
where the corresponding integrands for Jn, J12 and J13 are 
{i;+tiX,)-Ci{D,R)-^>{D,R)g,{i;+tTX^)-a{D:,R)-^{D,R)g,{i;+i^^^^ 

respectively. We use subindex b to denote the anomaly ('b' for bad). 

To estimate Jn, since the total signature is at most 3, thus the signature of ^{DyR)g 
is at most 1. Using Bianchi equations, up to an error of size CS*, we can replace the 
curvature term '^{D^R^g'^^ by V\l/g and we aslo replace the mild a{D^R) anomaly by a. 
Now we only consider the worst scenario. It should contain trx^. In this case, according 
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to signature considerations, V\l/g 7^ Va. Hence we control the worst scenario for Jn as 
follows, 

Jii < I [[ trx, • a ■ V*,| = I [[ trx„ • Va • < / ||Va|U2 (^^)||*,||^. 

JJV(u.u) JJV{u.u) Jo ^ ' ^ ' 

< n^{io + c(Xo)7^^ + cs's)< c(Xo)(7^ + + c^i 

To estimate J12, recall that the integrand is originally from (modulo the Hodge dual 
part) 

It has to be 

If ^ {(3,3),(3,c),(c,3),(c,d)|c,d = 1,2}, we know Di^I/ = 0. We also notice 

{/I, v) 7^ (c, d), otherwise Ji contains a double anomaly in curvature components. Hence, 
neither x or trx^ appears as connection coefficients. Furthermore, u is also absent. Thus, 
J12 is essentially K122 in last subsection (of the current section). Hence, 

Ji2<C5l 

To estimate J13, we need to carefully analyze the structure of 

DsRx'y' ■ D^L" ■ D.R^pzs, 

As we pointed out at the beginning of the section, the possible anomahes for last term 
are a{D4R), a{D^R), a{DsR) and /3{DcR). First of all, a{D4^R) can not occur since 
D^lJ" = 0. Thus, wc have three cases left to deal with. 

If '^{DyR) = a{DsR), thus u — 3 and the connection coefficient D^^l/ 7^ x or trx^. 
The integrand looks like 

Those terms were estimated in ^^122 in the last subsection. 

If '^{DyR) = a{D^R) or P{DcR), those anomalies are mild hence can be replaced 
by a. Combined with Bianchi equations if necessary, up to an error of size CS^, the 
corresponding integrand looks like 

ii^ + trXo) • V^'g • a. 

Those terms have been treated in Jn. Thus, we can derive estimates for J13. Combined 
with Jii and J12, we have 

Ji<c{Xo){n + nf^ + c6l (6.10) 
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6.3.3. Estimates for J2. In this case, both the curvature components in the integrand 
are anomalous. According to the different anomahes of D^Rx'^y^, we decompose J2 into 

J2 = J2I + ^22- 

For J21, D^Rx^Y^ = DzRiaAb] for J22, D^Rx^y^ = D^RsaSb- 

To estimate J21, we observe that DyR^pzs = D^RsaSb because the total signature is at 
most 3. Thus, up to an error of size C6^, 



J2I = I ^ ■ '^4" ■ Vsa = J2II + -^212 + -^213 

'V{u,u) 



V3UJ ■ Via ■ a+ u - a- V^a |^,~o + // u ■ a - VsV^a. 

V{u,u) J H^, JJv(u,u) 



For J2n, according to (2.5), up to an error term of size C64., we can replace Vso; by p, 
thus. 



We can then repeat the proof of Lemma |5.2| to derive 
Thus, 



||(p,«)llLf ^(H„)<c(Xo) + c(Xo)7^5+C53. 

{SC) 



J2ii<55 / (7^ + C(53)2(r5Xo + C(53) < +c(Xo) / n^du 
Jo Jo 

The last integral will be eventually removed by Gronwall's inequality, thus we could write 

J2u<C6l 



For J212, in view of Proposition 3.5, we have 



We derive 



II^IIl;* Ah^) < nn-^ + C6-^ < 7^3 + C6-^. (6.11) 

ysc) 

J Hu 

<55(7^! +C'53)(7^ + C'53)(r5c(Xo) + C53) < c(Xo)7^3 +C(53. 
For J213, we switch the order of V3 and V4 acting on a, 



J213 = I cu-a-V4V3a+ // w ■ a ■ [V3, V4]a 

'v{u,u) JJv(u,u) 
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We ignore the commutator term since it is bounded by CS^. The proof is routine. Thus, 
up to an error of size C6i, 

J213 = // V4C<; ■ a ■ Vaa + // 00 ■ V^a ■ V3Q; + / 00 ■ a - V^a \^~^ 

JJv{u,u) JJv{u,u) Jh^, 

The last two integrals are bounded hy CS^. In fact, we can replace the mild anomaly 
Vso; by a. Then the proof becomes routine. We now move on to the most dangerous 
term in J213, 

■ a ■ V^a 

It contains V4C0' which does not have L^^^^ estimates. We can only rely on its L^^^^ 

estimates. It requires an extra integration by parts. Up to an error of size C5^, we use 
(|2.10D to derive 

V4C0' ■ a ■ Vsa = // V4C0' ■ a ■ V/3 + // V^u-a-a 

V{u,u) jJv(u,u) JJv{u,u) 

The second integral is bounded by C5^, because we can put L^^^-j on V^uj and L^^^^ on 
curvature components. We move around the derivatives in the first integral, 

ViUj-a-V(3= [[ CO- V^a [[ u ■ a- V4V/3 + I w ■ a ■ V/3 

V{u,u) JJv(u,u) JJv{u,u) J H^i 

cu ■ a ■ V4V/3 + 

V{u,u) 

w-a- [V4,V]/3+ // w-a-VV4/3 

V{u,u) JJv{u,u) 



w-a- [V4,V]/3+ // Wuj-a-Wil5+ II w ■ Va ■ V4/3 

V{u,u) JJv{u,u) JJv(u,u) 

Thus, we removed all the anomalies. This yields the bound for J213. Together with J211 
and J212, we derive 

^21 < c(Xo)7^i + C(5i 
We now estimate J22- Its integrand is in the following form 

Since D^LT = Df'I^ = and D^R^aib is an anomaly, the integrand of J22 is reduced to 
the following form 

-Ds-Rs'^s* ■ -D'^L'^ ■ DcRdaib = -Ds-Rs'^a'' ■ x^^ ■ DcRdaAb 
We further decompose J22 into 

<^22 = <>^221 + <^222 
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such that the corresponding integrands of J221 and J222 are 



DsRs'^s' ■ ■ DcRdaAb-, D^Rs'^s" ■ trXg ■ DcRcaib, 

respectively. 

For J221, up to an error term of size C^J, we replace DcRdaib by V/3 + a to derive 



a b 



J, 



221 



X ■ Vsa ■ V/3 + // X ■ Vsa ■ a = J2211 + ^2212- 

V{u,u) JJv{u,u) 



For J2211, 

"^2211 



V{u,u) 



V{u,u) 



The last two integrals are bounded by C5i. The proof is routine. For the first term, we 

(-trx X - 2wx - a) ■ a ■ V/3 + Cb~^ < C6^ . 



use (2.2) to replace V3X, thus 

</2211 = 



For J2212, we have 



J'. 



2212 



V(u,u) 



Combining with(2.2) and (6.2), each of the above integrals generates the most difficult 
term 



X ■ a ■ a. 



In view of Proposition 3.5 



ML^^^iH) < c(Xo)r 4 + nini + C5-^ < c{io)5" + c{x^)n-^ + C5i. 



Thus, in view of Lemma 5.2 we have 

I / X-a-a| <^^(c(Xo)r^+c(Xo)7^t+C'53)(Xo + c(Xo)7^i+C(55)r2(Xo + 7^^+C5 

<c(Xo) + c(Xo)7^i+C(5l^. 
Thus, yields the estimates on J22i2- Combined with J2211, we have 

^221 < c(Xo) + c(Xo)7^^ + C5-^ 



I 16 



The last mission we need to complete is to show that J222 vanishes. In fact, the 
expression of the integrand 

-Ds-Rs^S^ ■ DcRcaAb = (^{D'iR)abSca l^{DcR)b 
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is not sufficient to show J222 vanishes. And we can not control it. In view of the formula 
(2.23), we have to also investigate the part from the Hodge dual to see the cancelation. 
We now show that 

Da^i^s'^s' ■ D*R,,^^ = a{D3*R)ai>eca*l3{D*R), = -aiDsR)abeca* PiD,R)b- 

We recall some basic formulas from [3]. For 1-form /3, *Pa = ^abP^- For a symmetric 
traceless 2-tensor a, *aab = £acC(^b and a*^ = aa^Ecb- One can easily check that*a = —a*. 
We also have =*a{W) = -a{Wy, /3(W) = -(3{Wy and Wab^c- Thus, 

a{D3*R)abeca*P{D:R)b = [-a{D^R)l^]e,aebdP{D*R)d = a{D:,R)l,ecnebd P{DcR)d 

= a{D^R)l,eebecaebdP{DcR)d = -a{D^R)l,e^5JP{D,R)d 
= -a{D:,R)abeca*l3{D,R),. 

Thus, J222 = 0. 

Remark 6.1. This cancelation can be explained in another way. To derive energy esti- 
mates at this situation, i.e. N = L and {X,Y,Z) = {L,L,Lb), by using Bel-Robinson 
tensors, it is equivalent to first commute (2. 16), (2. 14) and (2.17) with V3, and then mul- 
tiply the top order terms and integrate by parts directly on V{u,ub). In this way, we see 
directly that J222 type term do not appear. 

Putting all the estimates together, we derive 

J2<c(Xo) + c(Xo)7^3+C<5^. 

Combined with the estimates for Jq and Ji, we have 

J < c{Io) + c(Xo)(7^ + + 



Combined again with (6.1), (6.6) and (6.8), we derive 

+ II (V3P, Va^, V3/3) 11^2 ) < c(Xo) + c(Xo)(7^ + n)i 



||(V3/3,V3P,V3a) 



(6.12) 

Considering ||V3/3||i2 in (|6.12|), combined (|2.18|) and standard elliptic estimates for 
Hodge systems, we have 

II Va|U2 ) < c(Xo) + c(Xo)(7^ + + C5^. 

Considering V3P and V3a" in (6.12), combined (2.16), ( |2.14[ ) and elliptic estimates, we 
derive 



IIV^IL^ + IIV^IL^ (^J<c(Xo) + c(Xo)(7^ + ^)^+C5 

\SC.) — 



Hence, 



IIV^II^, .^(o.«). + ||(V^,V«)||^, .^(0,.) <c(Xo) + c(Xo)(7^ + E)« 

(sc) ^ > (sc) ^ — ^ ' 

We combine (Q, (Q, ( |5.12D and ( |6.13| ) to derive that 

7^ + ^<c(Xo)(l + 7^i) + CT^. 



(6.13) 
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In view of Section 2^, this completes the proof of the Main Estimates of the paper. 
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